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ABSTRACT 

Aims. The goal of this paper is twofold. Firstly, we present F- VIPGI, a new version of the VIMOS Interactive Pipeline and Graphical 
Interface (VIPGI) adapted to handle FORS2 spectroscopic data taken with the standard instrument configuration. Secondly, we inves- 
tigate the spectro-photometric properties of a sample of galaxies residing in distant X-ray selected galaxy clusters, the optical spectra 
of which were reduced with this new pipeline. 

Methods. We provide basic technical infomation about the innovations of the new software and refer the reader to the original VIPGI 
paper for a detailed description of the core functions and performances. As a demonstration of the capabilities of the new pipeline, 
we then show results obtained for 16 distant (0.65 < z < 1.25) X-ray luminous galaxy clusters selected within the XMM-Newton 
Distant Cluster Project. We performed a spectral indices analysis of the extracted optical spectra of their members, based on which 
we created a library of composite high signal-to-noise ratio spectra. We then compared the average spectra of the passive galaxies of 
our sample with those computed for the same class of objects that reside in the field at similar high redshift and in groups in the local 
Universe. Finally, We computed the "photometric" properties of our templates and compared them with those of the Coma Cluster 
galaxies, which we took as representative of the local cluster population. 

Results. We demonstrate the capabilities of F-VIPGI, whose strength is an increased efficiency and a simultaneous shortening of 
F0RS2 spectroscopic data reduction time by a factor of ~10 w.rt. the standard IRAF procedures. We then discuss the quality of 
the final stacked optical spectra and provide them in electronic form* as high-quality spectral templates, representative of passive 
and star-forming galaxies residing in distant galaxy clusters. By comparing the spectro-photometric properties of our templates with 
the local and distant galaxy population residing in different environments, we find that passive galaxies in clusters appear to be well 
evolved already at z ~ 0.8 and even more so than the field galaxies at similar redshift. Even though these findings would point toward 
a significant acceleration of galaxy evolution in densest environments, we cannot exclude the importance of the mass as the main 
evolutionary driving element either The latter effect may indeed be justified by the similarity of our composite passive spectrum with 
the luminous red galaxies template at intermediate redshift. 

Key words. Instrumentation: spectrographs - Methods: data analysis - Techniques: spectroscopic ~ galaxy clusters: general, redshifts 
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galaxy populations in clusters are dominated by red, passive 
galaxies, the most suitable feature to this aim is the continuum 
break at 4000A (D4000) and the associated absorption calcium 
lines (CaH/K). Forz ~ 0.8 - 1 this spectral region is redshifted 
to /I X! 7200 - 8000A, a wavelength window readily accessible 
for many ground-based telescopes. In particular, one of the cur- 
rently most efficient spectrographs able to cover this wavelength 
range and to simultaneously guarantee a very high and stable ef- 
ficiency of its CCD up to /I ~ 1 lOOOA is the FOcal Reducer and 
low dispersion Spectrograph (F0RS2, Sect. 13. Il l mounted on the 
UTl of VLT. Because of its performances at longer wavelengths, 
this instrument is widely used for imaging and spectroscopy of 
distant cluster galaxies and is hence the instrument of choice for 
the spectroscopic follow-up of the distant candidate systems in 
XDCP. 

At the time of writing XDCP provides the largest sample of 
X-ray selected distant galaxy clusters, with 30 confirmed clus- 
ters at z > 0.9 and a final aim of more than 50 clusters at 
z > 0.8 (30 at z > 1) to allow statistically meaningful evo- 
lution studies of the cluster population in at least three mass- 
and redshift bins. One of the main disadvantages connected to 
these expectations is the significant work related to reducing 
the large amount of data produced by the spectroscopic follow- 
up campaigns. Each XDCP target is observed with the F0RS2 
multi-object spectroscopy (MOS) mode, enabling an average of 
50 slits per mask. According to the survey expectations quoted 
above, the spectroscopic campaigns should finally yield several 
thousand slits, resulting in a similar amount of spectra to be re- 
duced and extracted (see Fig.|2] top). This process is very time- 
consuming if carried out with the traditional IRAF packages, and 
the pipelines provided by ESO for automating the procedures 
would not provide the necessary accuracy for an efficient extrac- 
tion of the spectra of such distant galaxies, which are mostly 
faint (I > 20I3. These reasons motivated us to develop a new 
dedicated pipeline for a quick and efficient reduction of spectro- 
scopic FORS2 data. Because of the many similarities between 
FORS2 and VIMOS data, we chose to build the new pipeline on 
VIPGI dScodeggio et al.ll2005l) and named it VIPGI. 

This paper is structured as follows: in Sect.|2]the main char- 
acteristics of VIPGI are briefly described to introduce the inno- 
vations of VIPGI, which are extensively discussed in Sect. [3] 
We then show an application of the new pipeline to a sample of 
distant clusters (Sect.|4]i that results in a new libray of spectro- 
scopic templates, while technical details and properties are dis- 
cussed in Sect.|5]and a conclusive summary is given in Sect.|6] 
Throughout this manuscript we assume a concordance ACDM 
cosmology, with Hq = 70 km s"' Mpc Qa = 0.7, Q.^ = 0.3 
and w = - 1 . 



2. The VIMOS Interactive Pipeline and Graphical 
Interface (VIPGI) 

VIPGI is a semi-automatic d ata reduction pipeline released in 
2005 dScodeggio et al.l l2005h written to process and archive 
the data obtained with the Visible MultiObject Spectrograph 
(VIMOS) mounted at the Melipal Unit Telescope (UT3) of the 
VLT in a quick and efficient fashion. It is able to handle the data 
taken with all the three available VIMOS modes: MOS, imaging, 
and integral field unit (IFU) spectroscopy. The creation of such a 
new, VIMOS-dedicated, reduction pipeline was motivated by the 



' Apparent Vega I-band magnitude of a z > 0.5 L* passively evolving 
galaxy with formation redshift z^„„„=5 and solar metallicity. 



need to process the huge amount of spectroscopic data expected 
from surveys such as the VIMOS -V LT D eep Survey (VVDS, 
Le Fevre et al. 2004) or zCOSMOS (iLiUv et al. 2007), each of 
which produced a total of ~50,000 spectra of galaxies in a very 
wide redshift range (0 < z < 5). 

The core of VIPGI is a library of C-written routines, th e 
VIMOS Data Reduction Software (DRS.IScod eggio et al.l2d0lh . 
currently used by ESO for the VIMOS online reduction. All 
these fundamental functions can interact with each other and 
with the user thanks to the adopted standard Python Tkinter 
graphical interface. This choice allows one to obtain a pipeline 
where the power and efficiency of the C code computation and 
the possibility of a continuous quality check by the user on the 
intermediate reduction results are both present. Namely, the user 
can constantly monitor the quality of the onging reduction step 
by step with VIPGI, repeating an intermediate reduction step if 
necessary without restarting the entire procedure. 

VIPGI was also designed to optimize the storage of a large 
amount of reduced data, with a clear and easily understandable 
filing strategy. More technical details on VIPGI are provided in 
the original paper of .Scodeggio et al., (,2005i) . 

3. The F0RS2-VIM0S Interactive Pipeline and 
Graphical Interface (F-VIPGI) 

Although VIPGI was specifically designed for the VIMOS in- 
strument, its capabilities are general enough to make it poten- 
tially usable with any other MOS spectrograph. This motivated 
us to consider VIPGI as a potentially useful tool for XDCP, 
whose cluster candidates have been (or are planned to be) spec- 
troscopically followed-up to safely confirm their nature of grav- 
itationally bound, distant systems. The instrument used for this 
purpose is the FOcal Reducer and low dispersion Spectrograph 
(FORS2) mounted on VLT. The XDCP spectroscopic campaigns 
have targeted ~70 cluster candidates since 2005, producing a to- 
tal of ~3,500 single spectra to be extracted and reduced. With the 
same original spirit of VIPGI we therefore adapted the power 
and efficiency of VIPGI to F0RS2 data as well. The result is 
F-VIPGI, a new pipeline that allows us to shorten the time for 
reducing FORS2 spectroscopic data by a factor of ~10 w.r.t. the 
standard IRAF procedures. 

3.1. The F0RS2 instrument 

FORS2 is the visual and near-UV FOcal Reducer and low- 
dispersion Spectrograph mou nted on the UTl unit (An tu) of the 
Very Large Telescope (VLT) (lAppenzeller et al 1 119981) . The in- 
strument covers the wavelength range from 330 nm to 1 100 nm 
with an image scale of 0.25"/pixel (or 0.125 "/pixel if the high- 
resolution collimator is used) in the standard readout mode (2x2 
binning) and a field of view (FoV) of 6.83'x6.83'. The detector 
consists of two MIT/LL CCID-20 chips, with 4096x2048 15yum 
pixels, each characterized by an excellent sensitivity toward the 
red part of the spectrum (up to A ~ 11 000 A) and the almost 
total absence of fringing pattern contamination. FORS2 can be 
used in many modes, including multi-object spectroscopy with 
exchangable masks, long-slit spectroscopy, imaging, spectro- 
polarimetry and high-time resolution imaging and spectroscopy. 

F-VIPGI was designed to work with all FORS2 data taken 
with the standard spectroscopic instrument configuration and 
straight slits, defined as those slits where the pixel-to-wavelength 
relation is "constant" along the slit length, and therefore the sky 
lines are perfectly aligned along the CCD rows (or columns, de- 
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Table 1. Properties of the grism forming the FORS2 standard instrument configuration for which F- VIPGI is usable. In parenthesis 
are the wavelength ranges actually used by F- VIPGI to provide the best-quality results. The listed values of resolution /1/A/l are 
computed at the central wavelength and for a 1 " slit. 



Grism name 


Central wavelength 


Wavelength range 


Dispersion 


Resolution 


Order separation 




[nm] 


[nm] 


[A/mm] 


[A/pixel] 


A/AA 


filter 


GRIS.600B+22 


465 


330-621 


50 


1.50 


780 


none 


GRIS.300V+10 


590 


330 (350) - 660 (925) 


112 


3.36 


440 


none 


GRIS.300V+10 


590 


445 (450) - 865 (850) 


112 


3.36 


440 


GG435-^81 


GRIS.300I+11 


860 


580 (600)- 1100 (1050) 


108 


3.24 


660 


none 


GRIS.300I+11 


860 


600- 1100(1050) 


108 


3.24 


660 


OG590+32 


GRIS.150I+27 


720 


330 (370) - 650 (980) 


230 


6.90 


260 


none 


GRIS.150I+27 


720 


445 (430) - 880 (990) 


230 


6.90 


260 


GG435+81 


GRIS.l 501+27 


720 


600 (590)- 1100 (1050) 


230 


6.90 


260 


OG590+32 



pending on the original orientation of the data). A summary of 
the properties of the FORS2 standard grisms equipment is pro- 
vided in Table[T] 



3.2. Conversion of F0RS2 files into fhe VIMOS format 

The functionalities of VIPGI were enabled not by creating 
new routines specifically coded for the FORS2 data format but 
instead by manipulating the FORS2 data themselves to convert 
them into the VIMOS format. Despite the differences between 
FORS2 and VIMOS CCDs architectures (with four squared sen- 
sors for the former and two rectangular ones for the latter), our 
idea was to make the software able to identify chips 1 and 2 
(hereafter Ql and Q2, respectively) of FORS2 as the first two in- 
dividual quadrants of VIMOS. In this way one can use the stan- 
dard VIPGI recipes with the "adapted" FORS2 fits files without 
affecting the quality and reliability of the data reduction results. 
The parts of FORS2 data that are manipulated by F- VIPGI are 
the image and the header. 



• image: the original FORS2 matrix of 4096 x 2048 pixels 
is transposed to obtain the same vertical orientation of the 
VIMOS frames. The wavelength dispersion orientation from 
lower (blue) to upper (red) side of the image is recovered 
in the same way. We also stress that this process does not 
affect the quality of the data because only the positions of 
the pixels are rearranged, whereas their information (counts) 
are left untouched. 

• header, the header structure of FORS2 files differs signif- 
icantly from those of VIMOS. In particular, while in each 
VIMOS file the header contains only the information on the 
slits of the chip to which the file belongs, in FORS2 each file 
has the information of all the slits of the mask on its header 
that are associated with Ql and Q2. During the conversion 
process F- VIPGI therefore changes these settings and re- 
moves all the information about the slits that do not belong 
to the chip of the frame itself. Only the six reference slits, 
common to both chips, are preserved along this process. We 
highlight that these changes are made for all FORS2 frames 
except the BIAS, as they do not contain any slit information 
in their headers. 



A specific and slightly different treatment is reserved to 
standard-star frames (STD), as discussed in Sect. 13. 71 



3.3. F-VIPGI calibration files 

As for VIPGI, the user has to initially provide specific calibra- 
tion files also for F-VIPGI. These are then used by the pipeline 
to locate the position and length of the slits on the CCD, to com- 
pute the coefficient of the interpolating polynomyal for the wave- 
length dispersion, to correct for bad CCD pixels, and to finally 
extract the single 1-D spectra. The calibration files are the fol- 
lowing: 

• grism table: it contains all the main spectroscopic infor- 
mation of the grism, i.e. central wavelength, wavelength 
range, and resolution. The emission sky lines usable for 
an additional refinement of the wavelength calibration (see 
Sect. l3.4l i are defined here, too. It is defined for each grism 
configuration. 

• CCD table: it maps all bad pixels/columns of the CCD that 
have to be corrected for by the pipeline. Two different files 
are needed for Ql and Q2, but these remain the same for any 
FORS2 observation. 

• line catalog: it contains the list of all arc lines used for 
the wavelength calibration (Sect. l3.4] i. For FORS2 only three 
files are needed, depending on the kind of lamp used (HeAr, 
HgCdHeAr or HgCdHeNeAr). 

• grism PAFfile: it is a parameter file in the standard ESO data 
interface control terminology. Basically, it is just a text file 
where information can be collected in a format that is very 
similar to the FITS keywords format and is then used to up- 
date FITS keywords in the header of FITS files. Specifically, 
this file contains all information related to the wavelength 
dispersion, optical distortion, and curvature of each CCD and 
each grism. 

• SPHOT table: it is a fits table containing magnitude vs. wave- 
length of the standard star that is used for the spectrophoto- 
metric calibration of the data. It has to be named exactly as 
the target name reported in the header keyword "ESO DBS 
TARG NAME". 

All the above calibration files except for the SPHOT table 
are provided together with the installation package and should 
not be modified by the user. 

3.4. Wavelength calibration 

One of the most powerful aspects of VIPGI, and hence of F- 
VIPGI, is that the user can perform sensitive calibration steps 
via visual tools that facilitate executing and constantly checking 
all the procedures. An example of this interactivity within the 
reduction process is shown in Fig.[T] Here the pipeline shows 
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Fig. 1. An example of a HeAr 
arc line catalog displayed on 
a raw lamp frame (in back- 
ground) within F- VIPGI. The 
green regions mark the ex- 
pected positions of the arc 
lines according to the infor- 
mation extracted from the fits 
header while the red ones on 
the right represent the line 
catalog used for the calibra- 
tion. The red part of the spec- 
trum is in the upper part of the 
image. 



the expected positions of the arc Unes of each sHt (green regions 
overlaid on the raw lamp frame) together with the line catalog 
used for the calibration (red regions on the right side). The user 
can iteratively refine the wavelength calibration across the entire 
chip by shifting each green line until the region line patterns 
perfectly match the underlying lamp frame. The positions of the 
line centroids are then computed automatically by the pipeline 
by applying a 2.5 sigma iterative rejection to remove lines that 
are deviate too far from the fit. 

This method produces a very accurate calibration with a 
median uncertainty on the wavelength calibration that roughly 
corresponds to one fifth of a pixel. As an example, for grism 
3001+ 11, whose spectra have a linear dispersion of 3. 24 A/pixel, 
the typical uncertainty on the wavelength calibration has an 

< rms >,i= (0.30 + 0.05)A, as shown in the bottom panel of 
Fig. 12] Another refinement of the calibration can also be obtained 
by using the position of some bright relatively isolated sky emis- 
sion lines whose wavelengths are known and defined in the grism 
table (see Sect. 13. 3b . For data taken with grism 3001+ 11 we used 
eight sky lines between A = 6300A and A = 10120A, which 
produced a final refined wavelength calibration with a typical 

< rms >,i < lA. This quantity can be converted into a corre- 
sponding redshift uncertainty of 6za ~ 1 - 2 • 10"'* by referring 
to Ac = 8600A. After the calibration process, all data are pre- 
liminarily reduced in the standard way by subtracting the bias 
frames, correcting for the CCD bad pixels and, finally, applying 
the flat-fielding correction. 



3.5. Sky line subtraction and atmospiieric absorption 
corrections 

A significant source of noise along the spectroscopic reduc- 
tion process is introduced by the strong absorption (telluric) and 
emission sky features produced by the O2 and H2O molecules 
and the OH" radical at ^ > 6000A (see Fig.[3] left). VIPGI 
is able to remove the two undesired sky contributions in a very 
efficient way, as shown in Fig.|4] The sky emission spectrum is 
computed for each slit in the free region.^ on the sides of the 
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Fig. 2. Top: Redshift distribution of 1543 FORS2 spectra pro- 
duced by XDCP and extracted with F- VIPGI. The solid line 
marks the secure assigned redshifts (>75% confidence on the 
assigned value) while the dotted line refers to unsafe ones (con- 
fidence <50%). The secure assessments represent ~12% of the 
entire sample at < z < 3 and increases to ~78% including also 
stars. Bottom: Distribution of rms wavelength calibrations com- 
puted for the corresponding 1230 reduced slits observed with the 
grism 3001+ 11. The median value is < rms >,)= (0.30 ± 0.05) A. 



science spectrum. If this area is large enough (at least 2" left 
between the spectrum and the slit edges, on both sides), the sky 
lines can be accurately modeled and efficiently subtracted from 
the final spectrum, as shown in the right panels of Fig.|3]and|4] 
After sky subtraction, we need to correct all the spectra also for 
the prominent sky absorption features (telluric lines) that could 



^ These regions are identified automatically by the pipeline with the 
following procedure: first a "slit profile" is traced by collapsing the data 



along the wavelength axis. The mean signal level and rms (sigma) of 
this slit profile are then computed with a r obust iterative proc edure by 
using the biweigth estimator described by iBeers et aP (Il990h . Finally, 
all groups of at least three pixels that are N sigma above the mean level 
are considered as object spectra while everything else is considered free 
sky region. In the above process the sigma threshold (A') can be set by 
the user. 
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8000 
Wavelength (A) 



Fig. 3. Left: Spectral absorption (top side) and emission (bottom side) features introduced by the night sky in the A > 6000A spectral 
region. Right: 2D emission spectra of the sky extracted for all slits of a FORS2 mask. The vertical alignment of the sky lines is 
indicative of a very good wavelength calibration over the entire mask. 
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Fig. 4. Example of the good quality achievable through the removal of absorption (left) and emission (right) sky features with 
F- VIPGI. Left panel: Atmospheric correction for a spectrum of a passive galaxy at z ~ 0.9 observed with grism 3001. Top: Flux- 
caUbrated spectrum (in black) without the atmospheric correction applied. It is evident that the telluric lines si A ~ 7600A strongly 
affect the S/N of the 4000-A break and the Call (H, K) lines. Bottom: The same spectrum after applying the atmospheric absorption 
correction. Now the Call lines are completely recovered. In red (solid and dotted curves) the positions of the most prominent telluric 
lines are shown. The black dashed lines mark the position of the most important spectral features in the displayed wavelength range. 
The best-fitting template of both spectra (an LRG one at z = 0.930) is overlaid in blue in both panels for reference. 
Right panel: Subtraction of the emission sky lines from two spectra with different redshift {spectmml at z = 1.086, spectrum! 
at z = 0.282). For each object the same part of the spectrum before and after the sky subtraction is shown at the top and bottom, 
respectively. In addition, the strongest spectroscopic features of each object are marked by the red boxes and are labeled on the top 
side. Although some residuals of the sky removal are visible in the bottom panels, the pipeline efficiently subtracts the undesired 
sky features and recovers the science features, which are sometimes completely outshone by the atmospheric emission. 



otherwise be interpreted as absorption lines in the object spec- 
tra. These spurious features can be removed provided that a suf- 
ficient number of well exposed flux calibrated spectra are avail- 
able. Even if the best results are obtained when the spectra se- 
lected for the correction span a wide redshift range (because the 
intrinsic features are efficiently removed in this way and only the 
sky ones are enhanced), a good telluric line absorption correction 
can also be gained by using cluster galaxies, whose redshift val- 
ues are relatively close to each other. If so, the pipeline finds the 



atmospheric coiTection by combining the spectra whose telluric 
lines are more evident and isolated and finally applies such a cor- 
rection to all slits of the mask. An example of the good quality 
achievable for telluric line absorption correction is shown in the 
left panel of Fig.|4l 
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Fig. 5. Image of the spectrum of a standard star taken with the 
FORS2 3001 instrument setup in the first chip (Ql). In red are 
shown the virtual slits of 30" length that F-VIPGI creates in the 
header of the scientific and calibration frames of the standard 
stars to mimic the VIMOS setting. In F0RS2 observations the 
standard star spectra are always imaged in Ql and are always 
enclosed in the first (#1) virtual slit. 




Fig. 6. Example of the spectrophotometric calibration for a FORS2 observation with grism 3001 and sorting filter OG590. Left 
panel: Extracted 1-D spectrum of the standard star, not yet corrected for the CCD response function. Right panel, top: The expected 
spectrum of the standard star provided by the SPHOT table (continuous line) and the fitted one (red points) according to the spectral 
resolution of the observation. The sensitivity function {bottom panel) is then computed as the ratio between the extracted 1-D 
spectrum and the fitted standard star spectrum converted into erg cm"^ s ' A ' . The user can edit the originally computed sensitivity 
function (red points) and choose the final one (continuous line) to apply to the single spectra for calibrating them into flux. 



3.6. The final product 

Once the preliminary reduction and the wavelength calibration 
are applied to each frame, these are calibrated in flux (see 
Sect. l3.7b and finally combined to obtain a stack of 2D sky- 
subtracted spectra. Single spectra are then extracte d from the 2D 
stacked frame, using a Horne optimal extraction (lHornelll986l) . 
and are saved in local folders in fits format. For each spectrum 
F-VIPGI also produces the associated noise vector generated by 
the extraction residual. 

3.7. Spectrophotometric calibration 

An important difference of the FORS2 observations to those of 
VIMOS concerns the methods used by the two instruments to 
acquire the spectra of spectro-photometric standard stars (STD). 
This step is mandatory to calibrate each spectrum in flux and 
to correct its continuum shape for the CCD sensitivity function. 
Standard stars in VIMOS are observed in all four quadrants and 
with a specifically designed mask composed of eight slits per 
chip, with fixed position and length (10")- In FORS2, instead, 
the STD frame is acquired only on Ql and by means of a single 
long-slit with 5" width. 

To overcome this discrepancy and enable the use of the stan- 
dard VIPGI routines for spectrophotometric calibration also for 
FORS2 data, F-VIPGI writes some keywords into the header 



of the standard star frames (and into the correlated calibration 
files) that mimic the presence of 11 slits in total (6 in Ql and 
5 in Q2j3with 30" length and 1" width. It is important to note 
that in FORS2 observations the standard star spectra are always 
centered on Ql at the position given by the header keywords 
CRPIXl and CRPIX2. The virtual slits are therefore created in 
the way to have the first slit (#1 of Fig.|5]l always centered on 
(CRPIXl, CRPIX2) and, so, always contain the star spectrum. An 
example of the part of the virtual mask relative to Ql is shown 
in Fig.|5] 

Once the standard-star spectrum is reduced and extracted, 
the pipeline computes the sensitivity function by comparing the 
observed stellar spectrum with the expected one contained in the 
SPHOT table (Sect. l3.3l l as shown in Fig.|6] The extracted curve 
is finally saved in a calibration table that can be applied to the 
1-D spectra for their flux calibration and for the correction of 
the CCD response function. Since standard stars in FORS2 are 
observed only on Ql, for Q2 data the same calibration table as 
computed for the first chip has to be used. An example of how 
the sensitivity function is computed for a standard star observed 
with the 3001 grism is shown in Fig.|6] 



' This is the maximum number of 30" length slits that can be uni- 
formly placed on the area of the two chips. 
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Fig. 7. Distiibution of the absolute values of the differences 
between the redshift measurements |Az| for 14 members of 
XMMUJ1230.3+1339 at z = 0.975 from multiple, independent, 
spectroscopic observations. The median value of < |Az| > = 
5.1 •10"'* is marked by the vertical dashed line. 

3.8. Redshift accuracy test 

The XDCP massive cluster XMMUJ1230 3+13 39 at z = 0.975 
dFassbender et al.l 1201 Ibt iLerchster et al.l 1201 ll) has been tar- 
geted by four different FORS2 spectroscopic observations in or- 
der to fully characterize its extended and dynamically unrelaxed 
galaxy population. This resulted in 65 confirmed spectroscopic 
members, 15 of which had multiple (double) observations. This 
enabled us to use these targets for a self-consistency test on the 
F- VIPGI redshift measurements, based on different independent 
spectroscopic observations of the same object. 

To this aim we selected 14 out of the 15 duplicated observa- 
tions, considering only those couples where the final extracted 
data were good enough to provide a safe redshift assessment for 
both cases. The median of the differences (in absolute value) 
between the redshift measurements is < |Az| >= 5.1 • 10^^, 
with a spread given by the semi-interquartile range (SIQ) of 
< |Az| >siQ^ 2.5 ■ lO-'*. This value can be translated into a cor- 
responding rest-frame velocity uncertainty of Sy^rest ~ 77 km s"' 
for a galaxy at z = 0.975. 

The distribution of |Az| is shown in the histogram of Fig.|7] 

4. An application to a sample of distant galaxy 
clusters 

In this section we show an application of f -VIPGI to the analy- 
sis of spectroscopic data for a sample of 16 distant XDCP galaxy 
clusters in the redshift range (0.65 < z < 1.25) observed with the 
grism 3001+1 1 (Table[T]). The main goal of this section is to pro- 
vide the community with a new set of spectroscopic templates 
that are intended to represent the galaxy population residing in 
distant galaxy clusters. For this reason, as shown in Sect. 15. II 
the present library is completely different from the existing ones 
(e.g., the ones constructed from SDSS, zCOSMOS or K20 sur- 



Fig. 8. Redshift distribution of the selected 187 galaxies, mem- 
bers of the 16 XDCP clusters reported in Table|2] 

veys), which contain templates representative of the population 
of field galaxies at different redshifts up to z ~ 1 .5. In addition to 
identifying redshifts of distant galaxies, the provided templates 
can also significantly improve the efficiency and the reliability of 
the photometric rtdshifl assessment o f distant clusters with c om- 
plex star formation histories (see, e.g.. lGuennou et al" and 
Pierini et a l. (2012) for a discussion) and they can be used to pre- 
dict rest-frame U - B colors of cluster galaxies at high redshifts. 

4. 1 . The spectroscopic sample 

For our study we selected all the XDCP clusters with at least 
three spectroscopically confirmed members with an S/N > 2 in 
the rest frame wavelength range of 4000 - 4300 A in their re- 
duced spectra. Each spectrum was also visually inspected to ver- 
ify that no contaminations from cosmic rays or bad sky subtrac- 
tion were present. For the cluster member selection, we firstly 
estimated the redshift cluster (zd) as the median of the redshift 
peak found in the proximity of the X-ray emission center After 
that, we selected the cluster members as those galaxies with a 
rest-frame velocity offset < 3000 km s"' from Zc/, corresponding 
to a redshift window cut of Az < 0.01 x (1 + Zd)- The list of 
the final targets selected for the next spectroscopic analysis is 
reported in Table|2] For each target a sequential ID, the spectro- 
scopic redshifts, the number of confirmed members used in the 
subsequent analysis, and the literature reference (if existing) are 
reported. A total of 187 cluster members have been identified in 
this way, with a redshift distribution in the range 0.65 < z < 1 .25 
as shown in Fig.|8] 

4.2. Spectrai indices analysis results 

We grouped the above galaxies into five diflFerent spectral classes 
according to the measurecQ values of the equivalent widths (EW) 

All spectral indices were measured in an automated way by means 
of a Python script developed by our group that uses the trapezoidal rule 
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Table 2. List of the 16 XDCP clusters used for our spectroscopic analysis. 



ID 


Zspec 




References 


cU 


0.7690 


13 




cl2 


1.2310 






cl3 


0.9410 


16 




cl4 


0.9590 




r70 in (?^ 


cl5 


0.8290 


8 




cl6 


0.7850 


5 




cl7 


0.6770 


12 




cl8 


0.8830 


4 


(1) 


cl9 


0.8940 


19 


(1) 


cllO 


0.7890 


13 


(1) 


clU 


1.1220 


7 


Cll in (2) 


cll2 


1.0740 


3 


(1) 


cll3 


0.8270 


9 


(1) 


cll4 


0.7460 


7 


(1) 


cll5 


0.9750 


52 


XMMU J1230.3+1339 in (3) and (4) 


cll6 


1.2030 


7 


(1) 



References. (1) To be published; (2) lFassbender et al.l ilOl lah : ('3^ lFassbender et al.l ll201 Ibh : W iLerchster et aLl(l2011l) 

Table 3. Classification criteria taken from lPoggianti et aD(l2009b used in this study for grouping the cluster galaxies into five pectral 
types according to the equivalent widths values of their [OII]/l3727 and H6 lines. The following convention is used: EW < for 
emission lines, EW > for absorption lines. The number of galaxies found in each class is given in the last column. 



Spectroscopic class EW([OII]) value EW(H<5/i) value Number 

[A] [A] of galaxies 

Passive >-5 <3 108 

Post-starburst >-5 >3 21 

Quiescent star-forming ]-25, -5] <4 23 

Dusty starburst ]-25, -5] >4 14 

Starburst <-25 any 21 



of their [Oil] and Hd lines, as summarized in Table|3] These 
two spectral features are, in fact, reliable indicators of ongoing 
and recent star formation activities within timescales of ~ 10^ 
and ~ 10^ yr for [Oil] and H5, respectively. The amplitude of 
their 4000-A break (D4000) and the equivalent widths of H/3, 
[OIII]/t4959 and [OIII]/15007 were also computed but were not 
used for the spectral classification. The errors on the indices 
were estimated by means of the noise vector provided by F- 
VIPGI for each extracted spectrum at the end of the reduc- 
tion process (see Sect. l3.6b . Specifically, for each spectrum 1000 
Monte-Carlo realizations were generated using the related noise 
vector as variance. The uncertainties on the indices were then 
computed as the rms of the results obtained by repeating the 
spectral indices analysis on the series of simulated spectra. 

We stress that the relative fraction of galaxies in the differ- 
ent spectral types is substantially biased by the method used in 
XDCP to se lect the targets for the spe ctroscopic follow-up. As 
discussed in iFassbender et al.l (1201 lah . to maximize the proba- 
bility of targeting actual cluster members, FORS2 masks are cre- 
ated in a fashion that the slits are preferentially placed on color- 
selected galaxies close to the expected red-sequence color and 
within the detected X-ray emission. This implies that the major- 
ity of the observed targets is actually expected to be passive and, 
hence, that no statistically significant conclusions about the rel- 



to perform the integration. The code functioning was tested on some 
spectra by comparing its results with those provided by the standard 
IRAF packages. 



ative number of galaxies in the different classes can be inferred 
from the last column of Table|3] 

We adopted the bandpasses defined in iBalogh et al.l (1 19991) 
for EW([OII]) and D4000 (hereafter D„40 00) while for EW(Hg) 
we us ed the definition of HSa given by IWorthev & Ottaviaiiil 
(Il997h . The galaxies were classified into five spectral classes fol - 
lowing the same method as described in lPoggianti et al.l (l2009l) . 
which is based only on the strength of [Oil] and H6a lines, as 
summarized in Table|3] 

Fig.|9] shows the distribution of the 187 galaxies in the 
EW([OII]) - EW(H5a) plane. The adopted loci for the five spec- 
tral classes are marked by the dashed lines and a color coding 
was used for a better visualization of the groups. 

The D„4000 index was then used to qualitatively test the va- 
lidity of our classification. The top panel of Fig. [10] shows the 
distribution of our galaxy sample in the EW(H5/i) vs D„4000 
space with the color code adopted in Fig.[TO]for identifying the 
different spectral classes. The distr ibution we find is remar kably 
similar to that reported in Fig. 1 of iGallazzi & Belli (l2009l) for a 
set of model galaxies. Namely, we observe that consistently with 
the results of these authors, the starburst (dusty and not) galaxies 
tend to be concentrated on the upper left side of the plot, indi- 
cating a relatively young (and still forming) stellar population. 
The passive and quiescent star-forming galaxies instead reside 
mostly in the center and on the lower right side, suggesting that 
they experienced the last episode of major star-formation activ- 
ity more than 3 Gyr earlier Finally, the post-starbursts are all 
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-5 5 10 

EW(H(5J [A] 

Fig. 9. Distribution of the 187 galaxies in EW([011]) - EW(H5a) 
plane. The dashed lines and the different colors mark the po- 
sitions of the fi ve spectral classes, with the same criteria as in 
iPoggianti et alJ (l2009l) . The color code is the following: blue: 
starburst galaxies; cyan: quiescent star-forming; magenta: dusty- 
starburst; red: passive; green: post-starburst. 



located above the sequence drawn by the passive and quiescent 
star-forming galaxies, a sign that their starburst activity ended 
within the previous 2 Gyr 

In the bottom panel of Fig.[TO]the distribution in the D„4000 
- EW([OIl]) diagram is shown. It is evident also here that for 
the adopted classification method the different spectral classes 
are confined to specific regions of the plots. In partic ular if we 
consider the relation adopted bv lFranzetti et al.l (l2007l) . 

D„4000 + EWi[Oll])/l5 > 0.7 

(marked by the dotted line in the plot) to divide passively evolv- 
ing (early-type) galaxies from star-forming (late-type) objects 
for a spectroscopic sample of V VDS galaxies at 0.45 < z < 1.2, 
all passive and post-starburst galaxies appear to be efficiently 
separated from the starburst ones. However, the plot also shows 
that the majority of quiescent star-forming and dusty-starburst 
galaxies are classified as early-type in this way. This effect was 
already predicted by the same authors, who claimed that some 
contamination in the early-type region was expected due to early 
spirals and that this is a general feature observed for color- or 
spectroscopic classification schemes. 

To identify possible active galactic nuclei (AGNs, e.g. 
Seyfert 2) erroneously ident ified as starburst galaxies, we used 

l2010l) . that was 



the diagnostic diagra m of iBongiorno et al.l 

originally defined in iLamareille et al.l (l2004h (and references 
therein), where the ratios between the EW of [OII]/l3727, 
HjS, and [OIII]/15007 lines are used to separate star-forming 
from AGN-photoionized spectra. The constraint of having the 
[OIII]/15007 line within the observed wavelength range 6000A 
< Aohs < 9800A makes this analysis feasible only for seven star- 
burst galaxies lying at redshift z < 0.95. 








^ -20 



S -40 - 




i- 



-60 



-80 



1.5 

D_4000 



Fig. 10. D„4000 versus EW([OII]) (top panel) and EW(H^a) 
(bottom panel) for the selected targets. The dotted line in the 
right panel represe nts the division betwee n early-type / late-type 
galaxies adopted in lFranzetti et al.l (l2007h . The color code is the 
same as in Fig.|9] 



In Fig.[TT] the distribution of these objects in the 
log(EW[OII]/l3727/EW(H;8)) vs log(EW[OIII]/l5007/EW(H;6)) 
plane is shown. Only one object (marked in red) appears to lie 
in the "AGN region" of the plot and was consequently excluded 
from the final stacking of the starburst spectra. 



5. A new library of spectroscopic templates 

Before proceeding with the final stacking procedure, we re- 
stricted the usable portion of each spectrum to the observed 
wavelength range of 6100 - 9280 A. This cut was made to avoid 
the higly sky-contaminated region at A > 9300A (Fig. [3] left) 
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Table 4. Spectral properties of the five averaged spectra of Fig. [13] The reported S/N values refer to the wavelength range 
4000A<i<4300A. 



Spectroscopic class 


S/N 


EW([OII]) 
[A] 


EW(H5/i) 
[A] 


D„4000 


EW(H^) 
[A] 


EW([OII1]/15007) 
[A] 


Passive 


9 


I.431±1.191 


0.148+0.981 


1.810+0.030 


1.829+0.606 


0.958+0.611 


Post-starburst 


6 


-0.178+1.739 


4.910+1.395 


1.543+0.034 


2.753+1.164 


0.540+1.448 


Quiescent star-foiming 


7 


-7.946±1.791 


0.593+1.231 


1.550+0.031 


1.664+0.887 


-1.184+0.985 


Dusty starburst 


7.5 


-7.719±1.605 


5.827+1.123 


1.412+0.030 


-2.357+1.849 


-5.842+0.693 


Starburst 


5 


-34.943±2.929 


1.368+1.887 


1.254+0.036 


-8.816+1.797 


-13.516+2.256 



o 
o 



^ 



1 
log(EW[OII]A3727/EW(H/3)) 



Fig. 11. Diagnostic diagram used bv lBongiomo et af] (1201 Ol) to 
distinguish between pure star-forming and AGN-contaminated 
galaxies applied to those galaxies of our sample classified as 
starburst (see Table|3) and lying at redshift 0.6 < z < 0.95 to have 
their [OIII]/I5007 line falling within the observed wavelength 
range. The solid and dashed curves show the demarcation and its 
+Q.15dex uncertainty, respectively, between pure star-forming 
galaxies (bottom region ) and AGN (top region), as defined in 
iLamareille et al.l (l2004h . The red point marks the starburst ob- 
ject possibly hosting an AGN that was therefore excluded from 
the final stacking procedure. 




4000 

Resl-frame wavelength [A] 



Fig. 12. Normalized rest-frame wavelength coverage for the 186 
galaxies used in the stacking procedure. The general result for all 
galaxies is represented by the thick solid curve, while the thinner 
lines indicate the different contributions from the five spectral 
classes defined in Sect. l4.2l The same color code as in Fig.|9]is 
used. 



while simultaneously preserving the observability of spectral 
features like H/3 and [OIII] for galaxies up to z ~ 0.85. 

The total rest-frame wavelength range covered by the com- 
posite spectra is 2700 - 5300 A, with a normalized wavelength 
coverage for all 186 considered galaxies shown in Fig. [12] In 
the same figure we also show the different contributions from 
the five spectral classes defined in Sect. 14. 21 We recall that here 
the possible AGN-contaminated starburst object of Fig. [TT] was 
excluded from our analysis. 

For the final stacking process each spectrum was rescaled to 
the mean flux computed in the wavelength range 4000A < A < 
4300A and then combined with the other objects of the same 
spectral class using the median as operator 

The resulting five composite spectra are shown in Fig. [13] 
For each spectrum the relative spectral class and the number of 
galaxies used in the stacking process are reported in the top left 



corner The most important spectral features recognizable in the 
rest-frame wavelength range AA - 2700 - 5300 A are labeled on 
the top side of the plot, and their expected positions are marked 
by dashed vertical lines. 

The spectral properties of the average spectra are quantified 
in Table[4| together with their S/N computed in the wavelength 
range 4000A <A< 4300A. 

We note that the values of EW(H5a) and, especially, D„4000 
obtained for the "passive" template are similar to those observed 
for the most massive and passive galaxies in the local universe 
by SDSS-DR4 ( p„4000|-=n ~ 1.9; EW (H^4)Un ~ -1.5) that were 
reported e.g. bv lGallazzi et al.l(l2005h . 

In the next two sections ( 15.11 and 15.2b we compare the 
spectro-photometric properties of our composite spectra with 
those of galaxies residing in diflferent redshift ranges and envi- 
ronments. 
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3000 3500 4000 4500 5000 

Rest-frame wavelength [A] 



Fig. 13. Final composite spectra, representative of the five considered spectral classes, created from a total of 186 galaxy spectra. 
The spectral classes and the number of stacked spectra are reported in the top left comer of each panel together with the expected 
position of the main spectral features marked by the blue vertical lines. A Gaussian smoothing filter of three pixels was applied. 



5. 1. Comparison of the resulting passive template with 
previous library spectra 

To highlight the original aspects of the spectroscopic library pro- 
vided in this paper and, hence, its uniqueness w.rt. other previ- 
ous similar works, we considered as an example the composite 
spectrum of our passive galaxies and qualitatively compared it 
with others representative of the same class of objects that reside 
in the local and distant universe and in different environments. 
Namely, we considered the distant sample of early-type galaxies 
in the field studied within the K20 s urvey, the composite spec- 
trum of which has been provided by iMignoli et al.l (l2005h . We 
highlight that the number of passive (stacked) spectra (93), their 



typical redshift values (< z > ~ 1.0, extended up to z ~ 1.25) 
and the spectro-photometric method used by the authors to iden- 
tify their passive galaxies are very close to the numbers and cri- 
teria adopted in our work (see Sect. 2). As a second reference 
sample we consider ed the nearby passive template provided by 
iKinnev et al.l ( Il996l) . produced from a total of four morpholog- 
ically selected bright (M^ <-21) elliptical galaxies residing in 
groups with < z >~ 0.0 08. Finally, we adopte d the compos- 
ite spectrum produced by lEisenstein et al.l (l2003h by using more 
than 1000 SDSS luminous red galaxies (LRG), selected in mag- 
nitude {Mg < -21.8) and in the redshift range 0.3 < z < 0.55, as a 
third element of comparison. 
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Fig. 14. Qualitative comparison of our composite passive spectrum (black line) with the empirical templates quoted in the text 
that is representative of passive galaxies residing in different environments and redshift ranges. All spectra are normalized in the 
wavelength range 3950A < A < 4050A and smoothed with a three pixel boxca r filter Left: The K2 composite spectrum of distant 
early-type galaxies in the field is shown in red while the template computed bv lKinnev et alJ ( 1 19961) for elliptical galaxies in groups 
at z ~ 0.008 is shown in magenta. It is evident that the passive galaxies of our sample are redder w.rt. the K20 ones at similar 
redshift, but bluer than the ellipticals in groups in the local universe. Ri^ ht: Comparison between the average spectrum of our 
passive galaxies and the LRG template provided bv lEisenstein et all (l2003h (in green) that was computed for z ~ 0.4, Mg < -21.8 
galaxies. Although the two spectra were obtained by using different galaxy selection criteria and redshift ranges, they appear to be 
remarkably similar overall (except for a slightly higher flux in U band for the LRG spectrum). 



All the above composite spectra, normalized in the wave- 
length range 3950A < A < 4050A and smoothed with a three 
pixel boxcar filter, are plotted in the two panels of Fig. [14] 

From this figure it appears that differences as well as similar- 
ities are present in our composite spectrum of passive galaxies 
residing in high-z dense environments and the other empirical 
templates considered here. In particular, we observe the follow- 
ing: 

1 . The passive galaxies in the densest environments at z ~ 1 ex- 
hibit a mean rest-frame U-V color that is redder than that 
of the coeval passive galaxies in K20, but bluer than that 
in groups at z ~ 0. Assuming that passive galaxies of the 
two distant samples span the same metallicity range, those in 
clusters seem to have evolved faster than their counterparts 
in the field. However, because passive galaxies in the XDCP 
are more luminous (massive) than those in K20, the differ- 
ence in their mean rest-frame U-V colors could be driven en- 
tirely by m ass and n ot (also ) by the density environment (cf. 
iThomas e t al. 20051 120T0I) . A more detailed spectroscopic 
analysis of our entire galaxy sample, aimed to characterize 
its star formation history, is in progress, however, and will 
be extensively discussed in a forthcoming paper 

2. The spectroscopic similarities between LRG and our passive 
galaxies would analogously suggest for the former sample 
an evolutionary path strongly driven by their mass and likely 
characterized by different assembly histories. 



In summary, we find significant differences of our passive com- 
posite spectrum compared to similar ones in dense low-z envi- 
ronments as well as high-z field environments. The closest re- 
semblance of our passive spectrum is found to the LRG tem- 
plate, obtained at intermediate redshift for very massive ellip- 
tical and bulge-dominated galaxies. All these findings seem to 
point toward a scenario where the evolutionary path of the pas- 
sive galaxies, also residing in clusters, is mainly driven by their 
mass rather than environment. 



5.2. U - B color of the templates 

In this section we compute the {U - B) rest-frame colors of the 
average spectra shown in Fig.[T3]and compare them to those typ- 
ical of nearby galaxy clusters. Specifically, we refer to the values 
observed for the Coma Cluster (z - 0.024) that were reported in 
the (U - B) vs H colo r-magnitude relation (CMR) of Fig. 3 in 
lEisenhardt eTal] (|2007|) . 

Firstly, we obtained the {U -B) rest-frame color [{U - B)iest] 
for each spectroscopic class by convolving the corresponding 
rest-frame spectrum with the response curves of the U and B 
filters of the KPNO 0.9m telescope. We also computed the 
quoted observed color by redshifting the average spectra to z = 
0.024 [(C/-B)-=o.o2] to obtain quantities comparable with Coma. 
These values are reported in the second and third column of 
Table|5] together with the differences between each rest-frame 
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Table 5. (U - B) colors of the five spectral templates at rest-frame and z - 0.024 (second and third column, respectively). For each 
class the differences between its (f/ - color and the corresponding value of a passive galaxy is reported in the fourth column. 



Spectroscopic Class 


iu - S)„,„ 


(U - fi);=0.02 


A((/ - S),„,| 


Passive 


0.553 


0.566 


0.000 


Post-starburst 


0.409 


0.455 


-0.144 


Quiescent star-forming 


0413 


0.432 


-0.140 


Dusty starburst 


0.229 


0.301 


-0.324 


Starburst 


-0.027 


0.041 


-0.580 



Table 6. The median apparent Vega H-magnitude (Col. 4) together with the number of used galaxies and their median redshift (Cols. 
2 and 3, respectively) for each spectral class. In Col. 5 the median of the difference between Hobs and the apparent H magnitude of 
a L* galaxy (at z = 0) passively evolving since Zform = 5 is reported. The errors refer to the semi-interquartile range. 



Spectroscopic Class 




< z > 


< ti„ts > 


< Hal,, -m*^> 


passive 


80 


0.94 + 0.07 


18.49 + 0.46 


-0.30 ± 0.47 


post-starburst 


11 


1.07 ± 0.12 


19.19 ±0.26 


-0.10 ±0.35 


quiescent star-forming 


14 


0.96 ± 0.02 


19.52 ± 0.39 


0.18 ±0.44 


dusty starburst 


18 


0.83 + 0.11 


18.71 +0.61 


0.14 ±0.39 


starburst 


10 


0.97 ± 0.14 


19.48 ± 0.55 


0.63 ± 0.89 




-0.5 - o - 

' o ° 

I 1 , 1 I 1 1 1 I 1 1 L I 1 a. 

-2 2 4 

H - m; 

Fig. 15. CMR of Coma Cluster galaxies (black pentagons) and 
the red sequence best fit, computed only for the early-type galax- 
ies defined as brighter than H - 14.5 (at the left of the dashed ver- 
tic al fine), with relative sc atter (black dotted lines) as provided 
bv lEisenhardt et all ( l2007h . The overplotted colored points mark 
the values obtained for the average spectra representative of the 
high-z cluster members, with the same color code as defined in 
Fig.|9l Note that for a consistent comparison all the observed Hi- 
band magnitudes were rescaled to the corresponding mj values, 
as explained in the text. For Coma galaxies a value of mj^^^^^ = 
11.13 was used. 

color and that typical of the passive class, taken as reference 
[A{U - B)resi\passive, fourth column]. 

After computing the five (t/-B);=o.o2 colors, we associated a 
median observed Vega H-magnitude (< Hoi,s >) to each spectro- 



scopic class to then plot these values in a color-magnitude dia- 
gram. We defined < Hohs > for each class as the median value of 
the apparent H-mag of the galaxies belonging to the class itself. 
However, only 13 out of the 16 XDCP clusters listed in Table|2] 
were observed in H-band, providing a total of 133 galaxies us- 
able to this aim. The number of the galaxies with H-band data 
for each class, their median redshift and Hoi,s are provided in 
Table|6l Finally, to compare the computed values with those of 
Coma galaxies in a consistent way, we rescaled all observed H 
magnitudes to a characteristic mj(z). It is defined as the appar- 
ent H magnitude of a galaxy at redshift z that has L = L* at z = 0, 
passively evolving since z/orm - 5. We adopted a suite of spec- 
tral energy distribution modelfl corresponding to the evolution 
of a simple s tellar population, computed withi n the framework of 
PEGASE 2 dFioc & Rocca-Volme range' 19971). The stellar mass- 
scale is matched to the observed mag nitude of an L* galaxy in 
clusters at z ~ in the Ks-band ( Str azzullo et a n i2006h . The as- 
sociation between age of the model and redshift of the observed 
galaxy is set by the adopted cosmological model. 

Since mj is a function of redshift, the appropriate value was 
computed and subtracted from Hots for each galaxy according 
to its redshift. The median < Hoi,s - > was finally ob- 
tained for each spectral class and is quoted in the fourth column 
of Table|6] The apparent H magnitudes of Coma galaxies were 
equivalently rescaled by referring to m HCnmn ~ 11-13 as quoted 
bv lEisenhardt et al.l (l2007h and found bv lde Propris et al.l(ll998h . 

The final (U - B);=o.02 vs Hohs - ni^ CMR for Coma and the 
average spectra of distant cluster members is shown in Fig. [15] 

From Fig.[T5]it is clear how the {U - B) colors of the Coma 
early-type galaxies are fully consistent with the typical {U - B) 
colors of the 0.6 < z < 1.2 cluster members classified as pas- 



In particular, these models assume a ISalpeteil (Il955h stellar initial 
mass function with lower and upper mass cut-offs equal to 0.1 and 120 
Mq, fixed solar metallicity and a range of 130 epochs with a time step 
of 0.1 Gyr since the instantaneous burst of formation. This choice pro- 
vides a suitable representation of the optical/near-infrared colors of ob- 
served, red, and dead galaxies associated with old, passively evolving 
stellar populations and an ellip tical morphology in different density en- 
vironments up to z ~ 2 (e.g., iPierini et aLll2004l : IWilman et al.ll2008l : 
iFassbender et al.ll201 lah . 
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sive, post-starburst, and quiescent star-forming from our analy- 
sis. This result may indicate that the stellar mass formation of 
the local passive cluster galaxies was almost completed already 
at z > 0.8, so that they did not experience any additional sub- 
sequent significant star formation process that could have made 
them bluer at z ~ 0. This appears to be confirmed because the 
values of EW(H5a) and D„4000 (reported in Tablegji of the com- 
posite spectrum of distant passive galaxies are similar to those 
found for the passive and most massive galaxies in the local uni- 
verse. These findings suggest that these galaxies have created the 
bulk of their stellar mass already at z > 0.8 and, according to the 
results discussed in Sect. l5.1l possibly even faster than the coeval 
galaxies in the field. As mentioned in Sect. l5.ll another study on 
the age-metallicity relation of the passive galaxies contained in 
our sample is currently underway and will be extensively dis- 
cussed in a forthcoming paper 

Finally, Fig. [15] also shows that the starburst galaxies at z < 
0.6 appear to be more luminous (in H-band) than the Coma 
galaxies with the same color. Even though a large part of this 
effect may be due to the expected selection bias toward the 
brightest objects in the distant universe, that in Coma very few 
star-forming objects appear to have H ~ in^coma ^ conse- 
quence of the downsizing scenario according to which the most 
vigorous star-formation activity at high re dshift was actually 
located in the most massive galaxies (e.g. , Cpwie et aL 19961: 



Gavazzi & Scodeggiol 19961 : iGavazzi et al.l [l996.: Juneau et alJ 
20051: IPannella et alj|2009r This effect is also beheved to have 



played a fundamental role in building-up the red sequence in 
galaxy clusters star ting from the more massive galaxies (see e.g., 
iRomeo et al.ll2008h . 



6. Conclusion 

We presented F-VIPGI, a new adapted version of VIPGI for re- 
ducing FORS2 spectroscopic data in a semi-automated and effi- 
cient way. We discussed the major improvements and technical 
aspects of the new software, which is now available to the scien- 
tific community, in the first part of the paper Additional infor- 
mation on the presented software are provided in Appendix A. 

In the second part we showed an application of the above 
pipeline to a sample of distant (0.65 < z < 1.25) X-ray se- 
lected galaxy clusters as part of the XDCP sample. We classi- 
fied all galaxies according to their spectral indices EW([OII]) 
and EW(H5) and stacked their spectra to create a new library of 
templates, which is particularly suited for a better redshift (spec- 
troscopic and photometric) measurements of galaxies residing in 
dense high-z environments. 

In the last part of the paper we compared the compos- 
ite passive spectrum with others representative of the passive 
galaxy populations residing in different environments and cos- 
mic epochs, finding some remarkable differences and similari- 
ties among them. 

Namely, we found that passive galaxies in clusters appear 
to be more evolved already at z ~ 0.8 w.rt. the field galaxies 
at similar redshift, supporting the idea that galaxy evolution is 
significantly accelerated in the densest environments. However, 
because our sample is significantly biased toward higher masses, 
we cannot exclude the importance of the mass as another driving 
element of evolution either We finally studied the (U-B) colors 
of the entire sample of composite spectra and compared them 
to those observed in the Coma Cluster The results of this work 
show that the colors of the passive galaxies in local clusters are 
fully consistent with those observed for distant cluster members 
classified as passive, post-starburst, and quiescent starforming 



according to our criteria. This suggests that the stellar mass for- 
mation in local cluster galaxies was almost completed already 
at z ~ 0.8. In addition, comparing the mean H-band luminosi- 
ties of the bluest galaxies of our sample with those in Coma, we 
found that at z > 0.65 the star-forming galaxies are much more 
luminous. Beyond the expected selection bias, the observed ef- 
fect may be partly explained by assuming a downsizing scenario, 
where the most vigorous starforming galaxies in the far universe 
were also the most massive ones. 



Appendix A: Additional information on F-VIPGI 

In the following appendix we provide some additional informa- 
tion useful to readers who consider using the software discussed 
in this paper. 



A.1. Compatibility witii tine current operating systems 

F-VIPGI is fully compatible with those architecture softwares 
where the old version of VIPGI can be installed and run. Since a 
few years, however, the plotting library Pmw.BLT is not updated 
anymore and hence started to become obsolete and incompatible 
with some recent operating systems and with all current 64-bit 
machines. 

In the following we list all operating systems that fully sup- 
port the current version of F-VIPGI with their version numbers 
reported in parenthesis: 

• Linux SuSE (9.x, lO.X, 11.0 and 11.1) 

• Linux Fedora (from 8 to 14) 

• Linux Ubuntu (from 8.xx to 10. xx) 

• Solaris (2.5) 

• Mac OSX (10.5 and 10.6) 

Another improved version of VIPGI is already being devel- 
oped and is expected to be completed by the mid of 2013. It 
will use a completely renewed graphical interface, fully usable 
also with the most recent architecture softwares, and it will also 
include reduction recipes suitable for LUCIFER data, the spec- 
trograph mounted at the Large Binocular Telescope. 

We finally stress that F-VIPGI is publicly released to the 
community as it is and without any guarantee of technical sup- 
port due to the lack of funds to support such assistance. 

A.2. How to obtain F-VIPGI 

The add-ons and the binary files of F-VIPGI can be downloaded 
from the following link: 

[littp://www.mpe.mpg.de/cosmology/cluster/FVIPGI| 

In the package, the cookbook of the new pipeline is also in- 
cluded. 



A.3. Link for spectroscopic templates 

The new library of spectra described in this manuscript, as well 
as some related technical information, are available in electronic 
form as ASCII and Flexible Image Transport System (fits) files 
at the CDS and also via 



http://www.mpe.mpg.de/cosmology/cluster/FVIPGI/specLibl 
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